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ABSTRACT 

We present cosmological parameter constraints from the SFI++ galaxy peculiar ve- 
locity survey, the largest galaxy peculiar velocity sa mple to date. The analysis is 
performed by using the gridding method developed in lAbate et al.1 (|2008l ). We con- 
centrate on constraining parameters which are affected by the clustering of matter: as 
and the growth index 7. Assuming a concordance ACDM model we find as = 0.9llo!i8 
and 7 = 0-55±g:i| after marginalising over f2 m . These constraints are consistent with, 
and have similar constraining power to, the same constraints from other current data 
sets which use different methods. Recently there have been several claims that the 
peculiar velocity measurements do not agree with ACDM. We find instead although a 
higher value of as and a lower value of VL m are preferred, the values are still consistent 
when compared with WMAP5. We note that although our analysis probes a variety of 
scales, the constraints will be dominated by the smaller scales, which have the small- 
est uncertainties. These results show that peculiar velocity analysis is a vital probe of 
cosmology, providing competitive constraints on parameters such as as- Its sensitiv- 
ity to the derivative of growth function, particularly down to redshift zero, means it 
can provide a vital low redshift anchor on the evolution of structure formation. The 
importance of utilising different probes with varying systematics is also an essential 
requirement for providing a consistency check on the best-fitting cosmological model. 

Key words: large-scale structure of universe - cosmological parameters - surveys - 
galaxies: kinematics and dynamics - galaxies: statistics 



J> 1 INTRODUCTION It is known that galaxies of different types cluster differently 

'kJI r™ , • , , • n , , ■ ■ n n , , • jDres slcr 1980: IConwav et a"fll2005l ). so they are clearly not 

1 the cosmological velocity field is induced by the gravita- , . , , . , , r , , . , , , . 

r '> , J J completely unbiased tracers ol the underlying mass, and this 

*H tional affect of the mhomogcncities m the matter distribu- . „ , , . , , . . „ 

rrt ■ ° issue is referred to gencrically as galaxy biasing. Compara- 

. . 1 tion in the universe. Therefore the recessional velocities of , . , ,, r .. , .. r , . 

tivefy therefore, the distribution of gafaxy peculiar veloci- 

galaxics experience deviations from pure Hubble now, the ,. , , ,. , ., ,. „ , . . 

. . . . „ ties has an advantage over the distribution of galaxies m 



result of the expansion of space, due to the distribution of 
matter both luminous and dark. The amplitude and distri- 



probing the matter distribution because they are likely to 

be unbiased traces of the matter velocity field. This is then 
bution of these deviations or peculiar velocities allows mea- . , , , , , , , , ., „ , , . n . „. 

simply related to the density held m linear theory, since pe- 
surement of the power spectrum of the fluctuations in the , ... , , r r,i 1 1 

cuhar velocities are a non-local function of the dark matter 
matter distribution. This is an important quantity to con- ,. , .. , . , n , , 

distribution then analysing the peculiar velocity held pro- 
strain because it tells us how structure has grown on differ- . , . „ , . , , , . . , 

vides information on scafes farger than the sampfed region 

ent scafes in the universe, and this m turn depends on the i TT „ ; — n l„„„-,k , ., . , . , . 

. , . . (Hoflman et al. 2001) as the velocity at a point is deter- 

amount of and nature of the components m the universe, i.e. . 7~, ~j . . , jn , . , 

mined by the mtegraf over the matter distribution m a farge 
the baryomc matter, dark matter, dark energy. , 

vofume. 

The distribution of galaxies in the universe is not likely 
to be the same as the distribution of matter, since most of Peculiar velocities are interesting for another reason as 

the mass is in the form of indirectly detectable dark matter. the y provide the only way to measure clustering of objects 

at practically redshift zero. Weak lensing, Lyman-a forest 
(Ly-a) and cluster count measurements are obtained only 
E-mail: abate@lal.in2p3.fr at higher redshifts, often where dark energy was just start- 
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ing to dominate. Understanding the strength and amount of 
clustering in the local universe could help discriminate be- 
tween different models that affect the growth of large scale 
structure such as dark energy models and modified gravity. 

In practice peculiar velocities are complicated by sev- 
eral factors. The major one is that on small scales the 
density field is highly nonlinear; these effects leak into 
the velocity field and cannot be described analytically. Re- 
sults fro m previous surveys which apply likelihood anal- 
ysis (see IZaroubi et ail Il997l . l200ll ; iFreudling et~aH 1 19991 . 
for Mark III, ENEAR and SFI respectively) seemed to 
overestimate the combination asO.^ 6 significantly compared 
to other probes at the ti me and to the current co ncor- 
dance cosmology . Stud ies bv lHoffman fe Zaroubil (|2000J ) and 
ISilberman et all <|200ll ) show that this over-estimation may 
be due to inaccurate modelling of the nonlinear part of the 
power spectrum, the small scales. The method employed in 
this paper overcomes this issue by smoothing the small-scale 
velocity field by averaging together peculiar velocity mea- 
surements regularly in a small volume of space, see Section 
13.21 for details. 

Another disadvantage of peculiar velocity analysis is 
the accuracy of the peculiar velocity measurements. Be- 
cause a redshift-independent estimate of the distance to 
the galaxy is required, and these estimates are usually 
extremely noisy. Most previous peculiar velocity surveys 
have used d istance indicators such as the Fundamenta l 
Plane/Dn-o- (jDiorgovski fc Dayisll 19871: iDressler et al.lll987h 
and the Tullv-Fisher relation l|Tullv fc Fisherlll977l '). These 
are galaxy scaling relations and it is because of their large 
scatter that the distance estimates are so noisy. The cal- 
culated distance is a relative distance measure which is 
strongly subject to a number of biases and also has very 
large uncertainties of around 20 percent, all of which tra ns- 
lates to the peculiar velocity fsee lStrauss fc W illick 1995|, for 
a review). Although type la supernovae (SNla) are more ac- 
curate as distance indicators (they have smaller distance er- 
rors, 5% to 10%) the current samples are very small (around 
150), whereas galaxy peculiar velocity samples are always 
in excess of 1000 obj ects. The SFI++ (|Masters et alj|2006t 
ISpringob~e t al. 2007) is the largest galaxy sample to date so 
provides an excellent opportunity to improve the statistics 
of peculiar velocity surveys. 

There has recently been a resurgence of inter- 
est in peculiar velocity analysis, after it fell out of 
favour in the late 1990s due to the apparent overesti- 
mation of os^m 6 described above. The renewed inter- 
est can be roughly catagorised into the following ar- 
eas: bulk flow measurements; as a "useful" systematic 
to SNla luminosity distance measurements; the kine- 
matic Sunyaev-Zeldo'vich (kSZ) ef fect. Bulk flow measure - 
ments have been recently mad e by lHaueb0lle et alj J2007h . 
IGordon. Land, fc Slosarl (l2008h.lFeldman fc Watkinsl (|2008f l 
and lWatkins. Feldman. fc Hudson! (2009). Measurements of 
the large scale velocity field, or bulk flow, has a number of 
interesting applications. If a difference exists between mea- 
surements of our dipole motion made using local peculiar 
velocities, and the dipole motion inferred from the cosmic 
microwave background (CMB) then this may be an indica- 
tion of a departure fro m the standard WMA P5 cosmological 
model (as published in lDunklev et al]|2009t) . Otherwise, be- 
cause the bulk flow is a sensitive probe of the matter density 



fluctuations on very large scales, it can probe the amplitude 
and shape of the matter power spectrum. 

Two recent studies o f peculiar velocit y surv eys in- 
dicated a large bulk flow. iKashlinskv et all (|2008l ) found 
a strong bulk flow on scales out to 300 /i _1 Mpc 
using peculiar velocity data fro m a cluster sample. 
IWatkins. Feldman. fc Hudson! |2009l ) reported a large-scale 
bulk flow beyond 50 /i _1 Mpc from a comprehensive compi- 
lation of peculiar velocity data. Both these papers claimed 
that the bulk flow at these depths, determined by their stud- 
ies, were difficult to explain within the framework of stan- 
dard ACDM model of cosmology. 

It has been recently found that the peculiar veloc- 
ity field does in fact affect measurements of the expansion 
of the universe from SNla. Low redshift (0.02< z <0.1) 
SNla are required both as a well-observed calibrating sam- 
ple and as a lever a r m for the SNla Hubb le diagram (see 
lAldering et~ail 120021 : lAbate fc Lahavl 120081 . for more de- 
tails). Although we are more interested in the expansion 
history of the universe at redshifts > 0.2, these low red- 
shift "calibrating" SNla have sufficiently low redshifts that 
the local velocity field systematically correlates their lumi- 
nosity distance m easurements. This effect must either be 
corrected f or (e.g. iNeill. Hudson, fc Conlevl l2007h or mod- 



elled (e.g. iHui fc Greene! l200rj : ICoorav fc Caldwelll 120061 : 



IGordon. Land fc Slosarl 12007) to retrieve an unbiased mea- 
surement of the dark energy equation of state w. 

The paper is organised as follows, in Section [2] we de- 
scribe the SFI++ data set, and in Section[3]we describe our 
method of analysing this data. Our results are presented in 
Section [4] and discussion and conclusions of these results are 
in sections [S] and HI] 



2 SFI++ PECULIAR VELOCITY DATA 

The SFI++ sample (Ma sters et"aLll2006T l consists of around 
5000 spiral galaxies that have measurements which are suit- 
able for the application of the I-band Tully-Fisher (TF) re- 
lation. This sample builds on the older SCI and SFI sam- 
ples, but includes significant amounts of new data, as well 
as improved methods for parameter determination. The ro- 
tation widths for the SFI++ sample come from the 21-cm 
line global profile widths (HI, about 60%) and optical ro- 
tation curves (Ha, ab out 40%). The catalogue presented in 
ISpringob et~atl (|2007f ) represents a homogeneously derived 
set of distances and peculiar velocities that are corrected for 
both the homogeneous and inhomogen eous Malmquist bi - 
ases using the 2MASS redshift survey l|Huchra et all 120051 ) 

for the line-of-sight density field. 

The data set we use comes from ISpringob et alj (I2009T ) 
after so me corrections were m ade to the original tables pub- 
lished in lSpringob et all (|2007l ). The morphological- type cor- 
rection for the peculiar velocity calculation was not properly 
applied in the original tables and this problem has now been 
rectified. The data set contains 4053 galaxies which have a 
median redshift of 0.0193. The magnitude of the error on the 
inferred distance is approximately 22% of the distance esti- 
mate. The sample covers most of the sky above the Galactic 



t The data set can be downloaded 
|http://www.iop.org/EJ/article/0067-0049/182/l/474/| 



from: 
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plane, with some deficiency of galaxies in the declination 
range of 5 = [-17.5°, -2.5°]. 



3 METHODOLOGY 

The a nalysis described in this Section follows lAbate et al.l 
(2008), for full details please refer to that paper. We use lin- 
ear theory to predict the velocity correlation function and 
use a multivariate Gaussian to calculate the likelihood. Sec- 
tion [372] describes how we overcome two potential problems: 
biases from nonlinear growth of structure and the large num- 
ber of velocities in the survey. 

3.1 Likelihood analysis using the VCF 

To estimate cosmological parameters we compute the radial 
peculiar velocity correlation function (hereafter VCF) from 
linear theory for each set of parameters. We define the pecu- 
liar velocity with observational error as Vi -$i = Vi = Si + 64, 
therefore the observed VCF is defined by 

Rij = {ViVj) = (SiSj) + (tiEj) (1) 

= + (2) 

where the average is over realisations of the universe. In 
practice we average over a large enough volume of space 
to assume ergodicity. The first term is the signal VCF and 
the second term is the contribution from the errors in the 
velocity measurements. Because the errors are assumed to 
be uncorrelated so they only affect the diagonal terms in 
Rij. In linear theory, the signal part ca n be split up 
into perpendicular and parallel com ponents l|Gorskil Il988l ; 
iGroth. Juszkiewicz. fc Ostrikerll989l ) , which are scalar func- 
tions of r=jr| 

£ij = cos0icos0j*n(r) +sin0iSin0j*x(r) (3) 

where the angles are defined by cos6*x = ? x • f and the 
diagonal elements £a are given by Eq. [5] below. The $n (r) 
and $i(r) represent the velocity correlations along the line- 
of-sight and perpendicular to the line-of-sight respectively. 
They are calculated from the matter power spectrum and 
assuming all galaxies are approximately at redshift zero 

*| U (r) = g ° 2/ 2 2 J" m) J P(k)B ll± (kr)dk (4) 

where B±_ = j (x)/x and Bn = j (x) and jo, jo are the first 
and second derivative of the zeroth order spherical Bessel 
functions respectively, Ho is the Hubble constant and fl m is 
the density of matter in the universe normalised by the crit- 
ical density, f(Q m ) is the derivative of the growth function. 
The auto correlation is given by 

bi = \ mf S m) Jnk)dk. (5) 

The dependence on the cosmological parameters of interest 
enters the above equations in the following ways: as through 
the normalisation of the power s pectrum P(k); fl m through 
f(fi m ) = W m (where 7 ~ 0.55, IWang fc SteinhardtJ Il998h 
and through its effect on the shape of the matter power 
spectrum; 7 through /(fi m ) as stated above. To compute 



Eq. [4] and Eq. [5]the power spectrum P(k) is generated using 
CAMB (|Lewis. Challinor. fc Lasenbvl[2000l ). 

The above equation for / assumes that the galaxies are 
at low redshift, therefore equations U and [S] are only valid for 
low redshifts. The full equation contains the growth rate at 
the redshift of each galaxy but because of the redshift range 
of SFI++ it is unnecessary to do the full computation here. 
Peculiar velocity surveys using distance indicator relations 
(as SFI++ docs) are unlikely to have a significant amount 
of data beyond a redshift of 0.05. The growth rate increases 
by less than 1 per cent between redshift zero and redshift 
0.05 for a flat ACDM model with fl m — 0.3, so this is a good 
approximation for this paper. We also use the approximation 
that the Hubble expansion is constant, described simply by 
a constant expansion rate for all galaxies in the survey. 

To calculate the covariance matrix R (Eq. [1] in matrix 
notation) it is necessary to calculate the equations[T]to[5]for 
a given set of cosmological parameters ©. Assuming that the 
peculiar velocities and the observational errors are Gaussian 
random fields the likelihood function for the parameter set 
can be written as 

L(&) = =L f cx P f -i V Vi (R _1 (©))ij VJ j .(6) 

7(27r) N |R(0)| V U / 

where Vi are the observed radial velocities. 

The likelihood analysis outlined above uses the indi- 
vidual galaxy peculiar velocities, Vi, as the data. However 
determining cosmological parameters in this way does not 
take account of the nonlinear part of the peculiar velocity 
signal because, as stated above, we make our prediction for 
the VCF based only on linear theory. The density field be- 
comes nonlinear only on small scales, above a wave number 
(fc) of about 0.2h Mpc -1 . The next section describes how 
we alter the above analysis to take account of this nonlinear 
signal. 



3.2 Gridding Method 

This method is a way of averaging together the peculiar ve- 
locities of spatially close galaxies by laying a grid across the 
survey. Averaging over a number of galaxies allows the linear 
signal to dominate. This averaging however will necessarily 
smooth the velocity field which is equivalent to damping the 
small scale contributions. This reduces the observed correla- 
tions because we average away some of the signal. If the data 
is averaged on a grid and inserted directly into the equa- 
tions in Section \3 . 1 1 without accounting for the averaging in 
the correlation function then the cosmological parameters 
will be biased. Therefore we present below how to imple- 
ment this type of binning on the velocity field, and detail a 
practical approach to take account of the binning accurately 
in the VCF. The technique we implement here is designed 
to be simple and fast. For a more complete description of 
this method, includin g its performance a fter rigorous testing 
with simulations, see lAbate et al.l (|2008l ). 

The method is implemented as follows. We lay down a 
grid across the survey and average together all the peculiar 
velocities within each grid cell so that 

v' m = (vi) iem (7) 



© 0000 RAS, MNRAS 000, 000-000 



4 Abate and Erdogdu 



{ah 



(8) 



where v' m is the radial peculiar velocity of the cell and e' m 
is the error on the velocity of the cell; the angle brackets 
denote an average over all galaxies i within the cell m. Note 
that €% is the contribution to the correlation function from 
the random velocity errors of each galaxy and therefore the 
remaining contribution e' m to the binned correlation function 
is reduced by the square root of the number of galaxies. 

This type of binning is then taken account of by mul- 
tiplying the power spectrum in Eq. [3] and [S] with a window 
function corresponding to the size and shape of the grid cell, 
which has the following form in Fourier space 



W{k) = 



8 sin [k x j) sin (fcy§) sin (k, 



k x 



k- 



(9) 



kefc 



where L is the length of a side of a cell in the grid, and the 
angle brackets denote an average over Fourier space direc- 
tions. This means equations that Eq. [4] and [5] become 



S mm 



Hp f_ ! (^ m ) 
2tt 2 

3 2tt 2 



W 2 (k)P(k)B ]u± (kr)dk (10) 



W 2 {k)P{k)dk 



(11) 



where we use the position of the cell centre to calculate all 
the required distances. The corresponding VCF and likeli- 
hood are then formed using the smoothed quantities. 

The method for accounting for the velocities described 
above assumes the data inside each grid cell is a continuous 
field, whereas it is in fact discrete values at the locations 
of the galaxies; we shall refer to this as the sampling effect. 
Galaxies trace discrete points of the peculiar velocity field, 
but if enough discrete points are averaged over then they 
will closely approximate averaging over a continuous distri- 
bution. However, some grid cells will not contain enough 
galaxies to provide a reasonable measure of the average of 
the velocity field within that cell, perhaps due to masked 
out areas in the survey or poor sampling in some areas. For 
example there is a deficiency of galaxies in the declination 
range of 8 — [—17.5°, —2.5°]. If this effect is unaccounted for 
it could bias the cosmological parameters, e.g. the effect on 
as would be to b ias it to too high a value. 

We found in lAbate et al.l (|2008T ) if we reduce the size of 
the modification to the diagonal elements of the VCF, £u, 
according to the number of galaxies in each cell we could 
remove any bias caused by the sampling effect. The following 
relation for this reduction was used: 



Ccorr 
mm 



+ 



(12) 



where is the corrected value used to calculate the di- 
agonal elements of the VCF; £,' mm is calculated from Eq. 1111 
£mm is calculated from Eq. [5] This correction uses the cor- 
rect value for the diagonal elements of the correlation func- 
tion in the limit that there is just one galaxy in the cell and 
also in the limit that there are infinite galaxies in the cell, 
i.e. a continuous field; for intermediate numbers of galax- 
ies the value of diagonal elements of the VCF should be an 
improvement on using 
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Figure 1. Best-fit <rg for various smoothing length L. The bold 
(blue) points are the results from the SFIH — h data, the faint 
(green) points are the results from the linear simulation. The 
dashed line is the fiducial as for the linear simulation. 



4 RESULTS 

We assume the following cosmological parameters through- 
out: h = 0.7, n s =l,Q b = 0.05 and fl k = 0. 

First we need to determine the optimum smoothing 
length L so as not to introduce bias in the cosmological pa- 
rameters from nonlinear signal or from the smoothing tech- 
nique itself. To do this we find the best-fit as for various 
smoothing length L for both the SFI++ data, and also a 
linear simulation of around 1500 galaxies. For this we have 
additionally fixed fl m = 0.3. The results of this are shown 
in Figure [T] The bold (blue) points are the results from the 
SFIH — h data, the faint (green) points are the results from 
the linear simulation. 

We simulate galaxy velocities for the linear simulation 
with the correlation function £ in Eq. [3] The number of 
galaxies in SFI++ is large enough that computation of the 
covariance matrix and eigenvalues was prohibitively slow 
so we simulate just 1500 galaxies. We use real positions of 
galaxies in SFI++ data set but simulate 50 realisations of 
peculiar velocities, and the resulting constraint on as is the 
average over these 50 realisations. The fiducial as for the 
linear simulation was as = 0.9 and is shown on Figure [T] by 
the horizontal dashed line. 

The large bias at L < 20Mpc//t is mainly an artifact 
of the gridding method, and not due in any major way to 
nonlinear bias or other systematic. This is confirmed because 
the linear simulation results follow closely the results of the 
SFI+- 1- data, and recover the fiducial as after L — 20Mpc//i. 
Therefore we find that L — 25Mpc//i is the best choice of 
smoothing length: large enough to avoid any bias but small 
enough so as not to redundantly increase the statistical error 
on the cosmological parameters. This value is close to the 
value of the optimum smoothing length found in lAbate et al.l 
(2008) of L = 20Mpc//i. 

With the optimum grid cell size of L = 25Mpc//z we 
constrain f2 m and as simultaneously. All cosmological pa- 
rameters except for fi m and as were kept fixed as specified 
at the start of this section, and the binning has been ac- 
counted for by using Eq. 1121 We vary fl m in two different 
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senses. First we vary the global value of Q m : the top panel 
of Figure [2] shows the 1 and 2-er likelihood contours in the 
0. m -as plane. Then we set Q m = 0.25 in the power spec- 
trum, and vary just the Q m in the derivative of the growth 
function f(Q, m ), which we shall now label Q f m . The 1 and 
2-cr likelihood contours of fl^-as are shown by the bottom 
panel of Figure [2] Both sets of contours have been overlaid 
on the joint constrai nts on erg and fl m fro m the 100 deg 2 
weak lensing survey (|Beniamin et al.ll2007h assuming a flat 
ACDM cos mology and adopting the nonlinear matter power 
spectrum of lSmith et al.l (2003). The weak lensing contours 
depict the 0.68, 0.95, and 0.99% confidence levels. The mod- 
els are marginalised, over h — 0.72 ± 0.08, shear calibration 
bias, and the uncertainty in the redshift distribution. 

The resulting constraints on fl m and erg are presented 
in rows 1, 2 and 3 of Table Q] Column 1 are the constraints 
when marginalising over (the global) Q m in the case of erg 
(as is marginalised over in the case of Sl m )- Column 2 is 
the same as column 1 except now these are the constraints 
where we have just varied Column 3 is where we set 
the global Q m = 0.25 (erg = 0.9 in the case of Q m ). Table 
[1] shows that no improvement on the constraints is gained 
when the effect of Q m in the power spectrum is fixed. All 
parameter constraints are consistent with ACDM and with 
WMAP5 to within 1-er confidence. 

At the end of Section 13.11 we mentioned the standard 
parameterisation of the linear growth function: 

/ S £i = nr (13) 
dm a 

where the growth index 7 = 6/11 ~ 0.55 for ACDM, 
and 7 ~ 0.55 for dark energy models with a s lowly 
varying equation of state dWang fc Steinhardtl Il998h . In 
the context of modified gravity the growth index param- 
eter can vary by as much as 30%, with the prediction 
for P GP branew orld gravity dDyali. Gabadadze. fc Porratl 
120001 ): 7 = 0.69 (|Linder fc Cahnll2007l ). Therefore observa- 
tional constraints on the growth index with better than 30% 
uncertainty begin to have the power to discriminate between 
dark energy and modified gravity models. 

In this spirit we also present the likelihood contours in 
the f2 m -7 plane in Figure [3] to see how the SFI++ peculiar 
velocity data performs with this goal in mind. Note that we 
only vary the global £l m simultaneously with 7. The contours 
shown are the 1 and 2-er, the dark/blue solid line shows 
the ACDM growth index (7 = 0.55) and the light/green 
dotted line shows the DGP growth index (7 = 0.69). After 
marginalising over fi m , the 68 percent confidence constraint 
on the growth index is 7 = 0.551^4 which is consistent at 
l-cr with Einstein gravity, and just consistent at l-cr with 
DGP gravity; see also the 4th row in Table [T] There is no 
constraint in column 2 because we do not vary 7 separately 
with Q.L.. 



Table 1. Parameter constraints from the SFH — h peculiar veloc- 
ities. Column 1 presents the l-cr constraint after marginalising 
over Q m (erg in the case of O m ). Column 2 is as column 1 except 
we have just varied H m in f(Q m ). Column 3 presents the 1-a 
constraint after setting Q m = 0.25 (erg = 0.9 in the case of Q m ). 



Parameter Margl Marg2 f2 m = 0.25 (or erg = 0.9) 




0.1 0.2 0.3 0.4 0.5 0.6 



£2 m 

Figure 2. The thick contours in both panels are the 1 and 2cr 
contours for the gridding method with a bin size of 25/i —1 Mpc. 
The top panel shows the constraints when the global Q m is varied, 
and the bottom panel shows the constraints when just the Q m in 
the derivative of the growth function is varied. The contours have 
been overlaid on the joint constraints on erg and Q m from the 100 
deg 2 weak lensing survey (grey scale contours) assuming a flat 
ACDM cosmology, depicting the 0.68, 0.95, and 0.99% confidence 
levels. 



5 DISCUSSION 

In this section we discuss our results in light of other recent 
constraints found in the literature. 
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5.1 Comparison with other constraints 

We jointly constrained erg and Q m , shown in Figure [2] where 
we have also plotted the constraints from the 100 deg 2 weak 
lensing survey. As can be seen from this figure, the weak 
lensing contours follow lines of constant crgf^, the normali- 
sation of the velocity signal. The contours from the peculiar 
velocities in the top panel however (where we vary the global 
fl m ) are more independent of the value of Q m . The joint ef- 
fects of Q m in the normalisation and in the shape of the 
matter power spectrum partially cancel out to produce the 
flatter contours. 

The bottom panel of Figure [2] shows the erg and fi£, 
1 and 2-cr contours, when just the fl m in /(f2 m ), i.e. the 
normalisation of the velocity signal divided by erg, is varied. 
The shape of these contours now matches the shape of the 
weak lensing contours; it follows lines of constant crgf^. 
This panel clearly illustrates how the different cosmological 
effects alter the shape of the contours. 

After marginalising over (the global) £l m we found 
erg = 0.91 ± 0.20, see Table [1] This is consistent wit h the 
value of erg estimated by Gordon. L and & Slosar ( 2007;) from 
SNla data. iGordon. Land fc Slosarl <|2007h used the VCF to 
calculate the peculiar velocity field induced covariance of the 
luminosity distances and, given a set of cosmological param- 
eters, fit this to the SNla data. They found erg = 0.79 ±0.22 
from a data set of 124 SNla with a median redshift of 
z — 0.024, similar to the median redshift of SFlH — h. It should 
be noted that their constraint includes priors on f2(,/i 2 and 
h, whereas we fix both fit and h to concordance values. 

These results are seemingly in confl ict with those found 
in lWatkins. Feldman. fc Hudson! ([2009) , who measured the 
bulk flow from several peculiar velocity surveys combined 
and found a bulk flow on large scales way in excess of what 
would be expected from ACDM. Their combined sample, 
while it is a combination of many different peculiar velocity 
surveys, is dominated by SFI+- h The bulk flow they mea- 
sured from SFI++ data alone was also inconsistent with 
ACDM. In addition they computed confidence levels in the 
crg-fi m /i 2 plane for the observed bulk flow (shown in Figure 
6 of their paper). They found that the WMAP 2-cr confi- 
dence levels for these parameters were excluded at greater 
than 2-cr confidence. Fixing all parameters except erg to cen- 
tral WMAP5 values they find erg ~ 1.7, with 2-cr and 3-cr 
lower limits of 1.109 and 0.878 respectively, excluding the 
WMAP5 value of erg at more than 3-cr. This result is a re- 
markably different to our constraint on erg after we fix all 
the other parameters to the concordance values stated at 
the start of Section U (these are only marginally different 
from the central WMAP5 values). We find erg = 0.9llo; 2 g 
and although the peak of the likelihood is at a relatively 
high value of erg, the result is still consistent with WMAP5. 

However our constraints are derived from a range of 
scales (the bulk flow is just the first moment of the flow on 
one scale) and these scales are mixed with one another, all 
contributing with different weights to the result. It is likely 
that in our measurement more weight is given to the 25-50 
h~ Mpc scale range because the larger scales are drowned 
out due to the larger uncertainties at greater distances. The 
measurement made bv lWatkins. Feldman. fc Hudson! (120091 ) 
was based just on the bulk flow on a 50/i _1 Mpc scale. The 
inconsistencies between our results and theirs may just be 



due to the different sensitivities of the bulk flow statistic and 
the VCF. 

As a check we perform a simple bootstrap test: we split 
the SFlH — h data into two equally sized samples, one sample 
contains the galaxies with the nearest measured distances 
and the other the furthest. We then repeat the analysis (con- 
straining the global Q m and erg simultaneously) on each sam- 
ple separately. We find no statistically significant difference 
between the two samples, and the marginalised constraints 



on erg are as follows: erg = 1.11_ 



(near sample), erg 



0.761 

0.27 (f ar sample). In fact we find the central erg value 
for the far sample, which probes mainly the larger scales, to 
be closer to the WM AP5 central value than the near s ample : 
an opposite result to IWatkins. Feldman. fc Hudson! f2009), 
but given the size of the uncertainties, this is not significant. 



5.2 Comparison with projected 6dFGS 
performance 

In lAbate et all (120081 ') forecasts were made for the perfor- 
mance of the S ix Degree Field Galax y Survey (6dFGS , 
iJones et al ]|2004T l galaxy velocity survev. lAbate et al.l (|2008T ; 
performed equivalent analysis to that done here on mock 
catalogues of the 6dFGS survey. They found the expected 
error on erg after marginalising over Q. m was approximately 
16 percent, compared to our current estimate of 20 percent 
from SFI++ data. This shows that the statistical improve- 
ment in this instance does not scale with the square root of 
the number of galaxies in the survey as might be expected. 
This can be understood by looking at the redshift distribu- 
tions of the 6dfGS mock catalogue and the SFI++ data: the 
median redshift of the 6dFGS mock is z — 0.040, a great 
deal larger than the SFI++ median redshift of z — 0.019. 
With a larger median redshift 6dFGS will have larger sta- 
tistical uncertainties on the peculiar velocities meaning less 
constraining power per galaxy than SFI++. 

The actual statistical improvement on parameter con- 
straints from 6dFGS over SFH — h will depend on the exact 
properties of the final data set. The release of the 6dFGS 
peculiar velocities is expected shortly (Magoulas et al. and 
Springob et al, in prep). 6dFGS wil l perhaps be most i mpor - 
tant for investigating the results of iKashlinskv et all (2008) 



and lWatkins. Feldman. fc Hudson! |2009l) . Because it 
pies such a large volume it will be possible to measure 
th e bulk flow at scales between 50/t~ 1 Mpc (as measured 
bv lWatkins. Feldman. fc Hudsonll2009f) and 300/i _1 Mpc (as 
measured bv lKashlinskv et al.ll2008l ) thus improving our un- 
derstanding of the bulk flow and confirming if there is a 
discrepancy with ACDM. 

5.3 Modified gravity constraints 

We also presented joint constraints on the growth index 
7 and Q, m , and found 7 = 0.55 ± 0.14 after marginal- 
ising over Q m , see Figure [3] and Table [1] This value is 
consistent with Einstein gravity and just consistent with 
DGP gravity and 1-er confidence. Other r ecent observa- 
tional cons traints of the growth index include Rapetti et al.l 



|2008l ) and[N csseris fc Perivolaropoulosl (|2008f ): we find our 
constraint to be consistent wit h these, and with a simi- 
lar amount of statistical error. iRapetti et all (120081 ) con- 
strained 7 from the X-ray luminosity function and found 7 = 
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°- 51 ±ai5 for tne ACDM model. Although they include ex- 
tra data sets (CMB, SNla, X-ray cluster gas- mass fractions) 
they also marginalise over a much larger parameter space, 
which includes allowances for s ystematic uncertainties, and 
so is a more rigorous analysis. iNesseris fc Perivolaropoulosl 
(2008) on the other hand constrain the growth function / in- 
directly, mainly through observations of ers(z) inferred from 
Ly-a, and they find 7 = 0.671°;?°. 

As the amount of data available for this type of analysis 
increases it will be interesting to see if the growth index re- 
veals any possible departures from the ACDM/dark energy 
predicted value of 7 ~ 0.55. 



6 CONCLUSIONS 

We presented cosmological parameter constraints using 
the SFlH — h galaxy peculiar velocity data. We performed 
the analysis by u sing the gridding method developed in 
lAbate et al.l (|2008ft . The gridding method is an efficient way 
to analyse the data given the large number of velocities in 
the survey, and without gaining biases on the cosmologi- 
cal parameters from the nonlinear growth of structure. By 
comparing with a linear simulation we found the optimum 
scale on which to average the peculiar velocity data was 
L — 25Mpc//i, and this scale is consi stent with the opti- 
mum scale found in I Abate et alj l|2008h . 
To summarise the main points: 

• Galaxy and SNla peculiar velocities provide the only 
way to truly measure clustering at redshift zero. They can 
provide a low redshift anchor for measuring the evolution of 
the growth of structure. 

• Peculiar velocities are sensitive to the derivative of the 
growth function, and can provide competitive constraints 
on parameters such as og and 7. Therefore they are a useful 
tool for discriminating between ACDM and modified gravity 
models. 

• Our constraints from SFI++ peculiar velocities are con- 
sistent with many other constraints found in the literature 
and with ACDM. However the constraints come from a 
mixed range of scales and necessarily will be dominated by 
the smaller scales which have the least uncertainties. Any 
departure from ACDM on larger scales (> 100/i _1 Mpc) will 
not be well tested by this data set. 

• The different degeneracy direction of constraints in the 
08-fi m plane compared to e.g. weak lensing will also be use- 
ful in breaking degeneracies, see also lGordon. Land fc Slosarl 
j2007ft . 

• Because peculiar velocities will be subject to different 
systematic effects to other probes such as cluster-counts, 
weak lensing or redshift distortions, they are an important 
complementary probe of cosmology. 

In the near future much more peculiar velocity data will 
become available, boosting even more the importance of 
this kind of analysis. The new technique of measuring 
galaxy cluster peculiar velocities using the kSZ could soon 
be achieved in large samples by Sunyaev-Zeldo'vich (SZ) 
surveys such as th e Atacama Cosmology Telescope (ACT, 
Fowler et all 120071) and the South Pole Telescope (SPT, 
Ruhl et al.ll2004 ). The kSZ is the temperature shift of CMB 
photons caused by a moving galaxy cluster. It is proportional 




Figure 3. 1 and 2a likelihood contours in the f2 m -7 plane. The 
dark/blue horizontal solid line shows the ACDM growth index 
(7 = 0.55) and the light/green horizontal dotted line shows the 
DGP growth index (7 = 0.69). 



to the line-of-sight momentum of the cluster gas, so therefore 
the radial peculiar velocity of the cluster can be inferred. The 
huge advantage of this method over using galaxies or SNla 
is that the kSZ measurement is redshift independent. This 
means we can probe the velocity field out to much larger red- 
shifts, and with a greate r accuracy than with galaxy pecu - 
liar velocity surveys, see iBhattacharva fc Kosowskv (|2008f ) 
for forecasts on the possible constraint on w. Although mea- 
surement of t he kSZ signal is still a difficul t challenge ob- 
serva tionally l|Knox. Holder, fc Churchll2004l ; iDiaferio et alj 
l2005h . 

On a shorter timescale, an upcoming release of the Six 
Degree Field Galaxy Survey will contain around 10 000 
galaxy peculiar velocities, an order of magnitude larger than 
most peculiar velocity surveys to d ate. Furthermore there 
are e ver growing SNla samples (e.g. |jha. Riess. fc Kirshnerl 
12007ft . and curren t and upcoming surve ys such as Skymap- 
per jKeller et all 120071 ) and GAIA (|Belokurov fc Evans! 
12003ft should push the sample size of SNla well into the 
thousands. 
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